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O 1 Abstract 

q_. Electric charge correlations were studied for p+p, C+C, Si+Si and 

centrality selected Pb+Pb collisions at -^/snn = 17.2 GeV with the 
NA49 large acceptance detector at the CERN-SPS. In particular, long 



range pseudo-rapidity correlations of oppositely charged particles were 
measured using the Balance Function method. The width of the Bal- 
' ance Function decreases with increasing system size and centrality of 

the reactions. This decrease could be related to an increasing delay of 
, hadronization in central Pb+Pb collisions. 
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1 Introduction 

Collisions of heavy ions have been used throughout the last decades in 
order to investigate the possible formation of the quark-gluon plasma (QGP) 
PQ, by studying a variety of characteristics 0. At the early stage of these 
collisions an extended region with large energy density may be produced, 
where hadronic may be replaced by quark-gluon degrees of freedom possibly 
leading to a new partonic phase of matter. In the subsequent evolution, 
the system dilutes and cools down, hadronizes and finally decays into free 
hadrons. These final state hadrons carry only indirect information about the 
early stage of the collision. 

Numerous observables, such as particle yields and measures of correla- 
tions and fluctuations, have been proposed that could signal the possible 
transition from the hadronic to the partonic phase. Recent data suggest 
that conditions, consistent with the occurrence of the QCD phase transition, 
are indeed reached [Hj in Pb+Pb collisions at 158^4 GeV at the CERN SPS. 
Moreover, results from the study of the energy dependence of single particle 
yields and spectra, suggest that a deconfined phase starts to be formed in 
the early stage of the reaction at low SPS energies jlj. The study of correla- 
tions and fluctuations is expected to provide additional information on the 
reaction mechanism of high energy nuclear collisions. In particular, event- 
by-event charge and mean pt fluctuations, have been analyzed over the past 
years Another important measure of correlations, the Balance Function 
(BF), was introduced by Bass, Danielewicz and Pratt [BJ. It measures the 
correlation of the oppositely charged particles produced during a heavy ion 
collision and its width can be related to the time of hadronization. The BF 
is derived from the charge correlation function that was used to study the 
hadronization of jets in p+p collisions at the ISR |2| and e~ + e + annihila- 
tions at PETRA [S]. The first results on the BF were obtained for Au+Au 
collisions by the STAR collaboration at RHIC 9J. 

In this paper we study the BF in p+p, C+C, Si+Si and centrality se- 
lected Pb+Pb collisions at a beam energy of 158^4 GeV, corresponding to a 
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center-of-mass energy of ^/saT/v = 17.2 GeV per nucleoli pair. The data were 
obtained with the NA49 detector at the CERN SPS. 

2 The Balance Function Method 

The motivation for studying the Balance Function comes from the idea 
that hadrons are produced locally as oppositely charged particle pairs. Par- 
ticles of such a pair are separated in rapidity due to the initial momentum 
difference and secondary interactions with other particles. 

Particles of a pair that was created earlier are separated further in rapid- 
ity because of the expected large initial momentum difference and the long 
lasting rescattering phase. On the other hand, oppositely charged particle 
pairs that were created later are correlated within a smaller interval Ay of the 
relative rapidity. Our aim is to measure the degree of this separation of the 
balancing charges and to find possible indications for delayed hadronization. 

In this paper the BF is used in order to examine the pseudo-rapidity (rj) 
correlation of charged particles. It is defined as a difference of the correlation 
function of oppositely charged particles and the correlation function of like- 
charge particles normalized to the total number of particles. The general 
definition of the BF reads j^j: 



where a and b could be different kinds of particles, whereas P\ and P 2 could 
be intervals in pseudo-rapidity. For example a could refer to all negative 
particles and b to all positive particles. Alternatively P 2 could be an interval 
of the relative pseudo-rapidity Arj = 1 77^, — 77^, | of the oppositely charged parti- 
cles, whereas Pi could be the interval of the pseudo-rapidity of the produced 
particles that is covered by the detector. In the numerator, N(b, P 2 \a, Pi) 
represents a conditional probability of observing a particle of type b in bin P2 
given the existence of a particle of type a in bin Pi. The terms N(b, Pi\a, Pi), 
N(a, P2\a, Pi), N(a, P 2 |6, Pi) and N(b, P 2 |6, Pi) are calculated using pairs 
from each event and the resulting values are summed over all events. For 
example, the term N(b, P 2 \a, Pi) is calculated by counting all possible com- 
binations of a positive particle in P 2 and a negative particle in Pi in an event 
and summing the number of combinations over all events. The other three 



P(P 2 |Pi) 



lr N(b,P 2 \a,Pi) -iV(a,P 2 |a,Pi) 
2^ N(a,Pi) 



N(a,P 2 \b,P 1 )-N(b,P 2 \b,P 1 ) 
N(b,Pi) 



(1) 
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terms are calculated analogously. The terms N(a,P\) and N(b,Pi) are the 
total number of negative and positive particles, respectively, that are within 
the studied pseudo-rapidity interval Pi, summed over all events. 

In our case, a and b are the negative and positive particles respectively 
that are within the pseudo-rapidity interval Pi and have a pseudo-rapidity 
difference Arj. So the definition of the BF takes the following form: 



= I v N + -(Ar,) - N-(A V ) N_ + (A V ) - N ++ (A V h 
{ 2^ AL N: )' { ] 



The most interesting property of the BF is its width. Early stage hadroniza- 
tion is expected to result in a broad BF, while late stage hadronization leads 
to a narrower distribution jHj . The width of the BF can be characterized by 
the weighted average (A77): 



k k 

(A V ) = Y,(B t -A Vt )/Y,B t , (3) 

8=0 1=0 

where i is the bin number of the BF histogram. 



3 Experimental Setup 

The NA49 detector JU] is a wide acceptance hadron spectrometer for 
the study of hadron production in collisions of hadrons or heavy ions at the 
CERN SPS. The main components are four large-volume Time Projection 
Chambers (TPCs) (Fig. [l]) which are capable of detecting 80% of some 
1500 charged particles created in a central Pb+Pb collision at 158^4 GeV. 
Two chambers, the Vertex TPCs (VTPC-1 and VTPC-2), are located in 
the magnetic field of two super- conducting dipole magnets (1.5 and 1.1 T, 
respectively), while the two others (MTPC-L and MTPC-R) are positioned 
downstream of the magnets symmetrically to the beam line. The set-up 
is supplemented by two Time of Flight (TOF) detector arrays and a set of 
calorimeters. The data presented in this paper are analyzed with a global 
tracking scheme jllj . which combines track segments that belong to the same 
physical particle but were detected in different TPCs. The NA49 TPCs allow 
precise measurements of particle momenta p with a resolution of cr(p)/p 2 = 
(0.3 - 7) • 10" 4 (GeV/c)- 1 . 
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The targets are C (561 mg/cm 2 ), Si (1170 mg/cm 2 ) disks and a Pb (224 
mg/cm 2 ) foil for ion collisions and a liquid hydrogen cylinder (length 20 cm) 
for hadron interactions. They are positioned about 80 cm upstream from 
VTPC-1. 

Pb beam particles are identified by means of their charge as seen by a 
Helium Gas-Cherenkov counter (S2') and proton beam particles by a 2 mm 
thick scintillator (S2). Both detectors are situated in front of the target. For 
p, C and Si beams, interactions in the target are selected by anti-coincidence 
of the incoming beam particle with a small scintillation counter (S4) placed 
on the beam line between the two vertex magnets. For p+p interactions 
at 158 GeV this counter selects a (trigger) cross section of 28.5 mb out of 
31.6 mb of the total inelastic cross section. For Pb-ion beams, an interaction 
trigger is provided by anti-coincidence with a Helium Gas-Cherenkov counter 
(S3) directly behind the target. The S3 counter is used to select minimum 
bias collisions by requiring a reduction of the Cherenkov signal by a factor of 
about 6. Since the Cherenkov signal is proportional to Z 2 , this requirement 
ensures that the Pb projectile has interacted with a minimal constraint on the 
type of interaction. This setup limits the triggers on non-target interactions 
to rare beam-gas collisions, the fraction of which proved to be small after 
cuts, even in the case of peripheral Pb+Pb collisions. 

The centrality of a collision is selected (on-line for central Pb+Pb, Si+Si 
and C+C and off-line for minimum bias Pb+Pb interactions) by a trigger 
using information from a downstream calorimeter (VCAL), which measures 
the energy E of the projectile spectator nucleons. 

4 Data Analysis 

4.1 Data Sets 

The data sets used in this analysis come from p+p, C+C, Si+Si and 
Pb+Pb collisions at 158^4 GeV. For Pb+Pb interactions data with both 
central (2 • 10 5 ) and minimum bias trigger (6 • 10 5 ) have been analyzed in 
order to study the centrality dependence of the BF. The minimum bias data 
were subdivided into six different centrality classes |T2] according to the 
energy recorded by the VCAL, from class Veto 1 (the most central collisions) 
to class Veto 6 (the most peripheral collisions). The most central Pb+Pb 
interactions correspond to 5% of the total geometric cross section (Table GJ. 
Since minimum bias data provide only a small number of central collisions, 
we used in addition trigger selected central data. Finally, we analyzed three 
different data sets (Table |2J) of Pb+Pb minimum bias events coming from 
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two different data-taking periods (1996 - Data Set 1 and Data Set 2, 2000 - 
Data Set 3) with opposite magnetic field polarities (positive field polarity - 
Data Set 2 and Data Set 3, negative field polarity - Data Set 1) in order to 
estimate the systematic uncertainties (Table EJ). 

The event centrality is characterized by the mean impact parameter (b) 
and the corresponding number of wounded nucleons (N\y)- For each bin 
of centrality these quantities were determined by use of the Glauber model 
as implemented in the VENUS event generator [13j . In order to estimate 
the correlation between the energy deposited in the VCAL and (b) or (N\y) 
minimum bias VENUS events were processed through the GEANT detector 
simulation code, and the energy deposited in the VCAL was simulated. All 
these quantities are listed in Table [TJ 

4.2 Event and Track Selection 

In order to reduce the contamination from non-target events and non- 
vertex tracks, selection criteria were imposed both at the event and the track 
level. 

Events were selected that had a proper position of the reconstructed pri- 
mary vertex. The vertex coordinate V z along the beam axis had to fulfill 
\V Z — V ZQ \ < Az where the values of the central position V ZQ and the range 
Az are shown in Table El for p+p, C+C, Si+Si and Pb+Pb reactions, re- 
spectively. In addition the vertex coordinates V x and V y perpendicular to 
the beam axis had to fulfill \V X — V xo \ < Ax and \V y — V yo \ < Ay, where the 
values V Xo ,V yo and Ax, Ay can also been seen in Table El for all the data 
samples analyzed. 

Selection criteria at the track level were imposed in order to reduce the 
contamination by tracks from weak decays, secondary interactions other 
sources of non vertex tracks. Thus, an accepted track had to have an extrap- 
olated distance of closest approach d x and d y of the particle at the vertex 
plane within the range: \d x \ < 2.0 cm and \d y \ < 1.0 cm. In addition the 
potential number of points in the detector for the selected tracks had to be 
more than 30. To suppress double counting due to track splitting, the ratio 
of the number of reconstructed points to the potential number of points was 
required to be larger than 0.5. 

The NA49 detectors provide large acceptance in momentum space; how- 
ever the acceptance in the azimuthal angle <fr is not complete. The boundary 
of the acceptance region can be described with the formula [Ti] : 

^)-^T(W + B - (4) 
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where the values of the parameters A, B, C and D depend on the rapidity 
interval and are given in Table El (see Fig. E] for examples) . The inclusive 
pseudo-rapidity distribution after applying the acceptance filter can be seen 
in Fig. El 

Finally we required tracks to additionally satisfy the following criteria: 
0.005 < p T < 1.5 GeV/c and 2.6 < 77 < 5.0. As shown in Fig. El the 
phase space analyzed covers most of the forward rapidity region, where the 
geometric acceptance is maximal. 

4.3 Results 

In this section, we will present results on the BF (Eq. EJ) measured in 
p+p, C+C, Si+Si and Pb+Pb at ^/s^n = 17.2 GeV that were subjected to 
the event and track quality as well as to the phase space cuts described in 
the previous section. 

In order to study the centrality dependence of the BF, we analyzed 
Pb+Pb collisions that were divided into six centrality (Veto) classes [T2"] . 
from 1 (the most central collisions) to 6 (the most peripheral ones) (Table 

®- 

The results are shown in Fig. 0] where the BF is plotted as a function 
of Arj, the pseudo-rapidity difference of the charged particles. The error 
on each measured point is the statistical error. For visual comparisons the 
distributions were fitted with a Gaussian function having a fixed mean at zero 
(curves in Fig. @J. From inspection of Fig. 01 as well as from the values of the 
weighted average (Arj) that are listed in Table IH we notice that the width 
(A77) of the BF is narrower for the most central collisions (Veto 1) than for 
the peripheral ones (Veto 6). It should be mentioned, that for the calculation 
of the width (Eq. EJ) we excluded the first point of each distribution, since 
from ^5] it was shown that this point is significantly influenced by Coulomb 
interactions and Bose-Einstein correlations. 

Furthermore, in order to extend the method to a system size study, we 
have analyzed C+C and Si+Si collisions at ^snn = 17.2 GeV. The BFs for 
the data samples of these two systems are shown in Fig. El The distributions 
are wider than those of the most central Pb+Pb collisions and tend to be 
similar to the ones coming from the most peripheral Pb+Pb interactions 
(Veto 6). This conclusion is confirmed by the corresponding values of (Arj) 
displayed in Table |H 

In addition, we have studied p+p interactions at y / s/v r /v = 17.2 GeV. The 
resulting BF distribution shown in Fig. El is significantly wider than that for 
Pb+Pb interactions. The calculated widths {Arj) for p+p, C+C, Si+Si and 
all centrality classes of Pb+Pb interactions are summarized in Table E] along 



8 



with their statistical errors. 

4.4 Systematic errors 

The systematic errors of the width of the BF were estimated by varying 
the cuts in V z , d x and d y and by comparing results obtained from different 
data taking periods. The results are described in this section. 

The dependence of the width of the BF on the cut Az for the event vertex 
position and the upper limit cuts on the impact parameters \d x \ and \d y \ are 
shown in Fig. El for p+p and Pb+Pb (central and peripheral) collisions. 
The resulting variations of the width of the BF are used to estimate the 
systematic errors due to contamination of non-target interactions and non- 
vertex tracks. They amount to no more than 0.006, 0.009 and 0.003 for p+p, 
Pb+Pb peripheral and Pb+Pb central collisions, respectively. 

Finally, as mentioned in a previous section, we analyzed three different 
data sets of minimum bias Pb+Pb collisions. The observed difference in the 
BF width are smaller than: 0.005, 0.009, 0.006 and 0.004 for Veto 6, Veto 5, 
Veto 4 and Veto 3 centrality selection respectively. 

To summarize, the estimated systematic errors of the width of the BF for 
p+p, C+C, Si+Si, Pb+Pb peripheral and Pb+Pb central collisions are no 
more than : ±0.006, ±0.010, ±0.012, ±0.009 and ±0.003, respectively. 

5 Discussion 

In this section the results presented in the previous ones, will be compared 
to models and to results from RHIC obtained by the STAR collaboration [H] . 

The BF for each centrality class was calculated for mixed events that were 
produced by randomly choosing particles from different events with similar 
vertex position and multiplicity. As shown in FiglH the BF for mixed events 
goes to zero because of the removal of correlations caused by global charge 
conservation. Another method of mixing was applied to the data sample in 
order to estimate the maximum possible value of the width of the BF while 
retaining the constraint of charge conservation. This shuffling procedure [§| 
is a mixing method in which the value of the pseudo-rapidity of each track 
is taken randomly from the collection of pseudo-rapidity values of the tracks 
in the same event, whilst keeping the charge of each track the same. The BF 
for shuffled data is broader for each centrality class than the one obtained 
from the real data (Fig. HJ). The values of (A77) for the shuffled data analysis 
are listed in Table El 
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Finally, in order to further investigate the origin of the system size and 
centrality dependence of the BF, we generated p+p, C+C and Si+Si collisions 
as well as centrality selected Pb+Pb interactions at y / J/v r /v = 17.2 GeV using 
the HIJING event generator PB]. The model is based on the excitation of 
strings and their subsequent hadronization according to the LUND model. 
The latter contains short range correlations of oppositely charged hadrons 
which are consistent with measurements from e + + e~ annihilations. The 
rescattering of produced hadrons is not included in the model. 

The generated data sets were analyzed with and without applying the 
NA49 acceptance filter. The results revealed that the acceptance filter slightly 
increases the width by about 4%. This suggest that this filter removes a frac- 
tion of balancing charges. The filtered distributions for Pb+Pb collisions and 
interactions of lighter systems are plotted in Fig. 0] and Fig. |3] respectively. 
The values of the widths are included in Table El The BF for HIJING is in- 
dependent of centrality and system size and is wider than the one calculated 
from the real data for central, mid-central and mid-peripheral collisions. On 
the other hand both HIJING and real data distributions tend to be similar 
for the most peripheral Pb+Pb collisions (Veto 6) as well as for the lighter 
systems. 

In order to demonstrate the dependence of the BF 's width (At]) on the 
centrality class in Pb+Pb interactions, the BFs in different centrality bins 
were normalized to the same area and plotted on the same graph (Fig. |7j). A 
significant narrowing of the BF width with increasing centrality is observed. 

Fig. |H1 shows the dependence of the width (A77) of the BF on the mean 
number of wounded nucleons (Nw) (Table HJ). The results for p+p, C+C 
and Si+Si collisions are also included. The width decreases monotonically 
with (Nw). On the other hand the width of the BF from both HIJING and 
shuffled data does not show any clear dependence on centrality. 

Fig. El shows the dependence of (A77) on the normalized mean impact 
parameter (b)/b max . The values of the impact parameter are listed in Table 
H Once again the strong decrease of the width with increasing centrality of 
the collision is obvious. The results from a similar analysis performed for 
Au+Au collisions at y / J/v r /v = 130 GeV by the STAR collaboration at RHIC 
jH] are plotted in Fig. The width of the BF decreases from peripheral 
to central collisions by 17 + 3% for the NA49 data, whereas for the higher 
energy STAR data the corresponding decrease is of the order of 14 ± 2%. 

The narrowing of the BF compared to shuffled events is of similar mag- 
nitude in both experiments. The somewhat smaller difference between the 
widths for data and shuffled events for NA49 may be due to the incomplete 
azimuthal acceptance. 
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The influence of the decay of resonances on the width of the BF was 
estimated using the HIJING event generator. We found that the BF width 
increases by about 4% when p°-meson decays are switched off. In the model 
the fraction of pions coming from p° decays (about 19%) is approximately 
independent of centrality. Therefore, the effect of p° decay can not explain 
the strong system size and centrality dependence of the width of the BF that 
we observe in our experimental data. 

The measured narrowing of the BF is qualitatively consistent with the 
delayed hadronization scenario jfH IHj of an initially deconfined phase. Several 
model calculations have been published which provide a more quantitative 
description ^3 EE EH| • In particular, within models based on statistical 
hadronization and hydrodynamic expansion the width of the BF was found to 
decrease with increasing transverse collective velocity of the matter at freeze- 
out [I3[TB1 EH] an d thus with the collision centrality. However, a quantitative 
description of the STAR data was possible only when the condition of global 
charge conservation (a single fireball model) [THJ EH] was substituted by a 
stronger condition of charge conservation in sub volumes (a multi-fireball 
model) |17j . The quark coalescence model was applied to the hadronization 
of the deconfined phase in [20J . When including radial flow, good agreement 
with the STAR measurements was obtained also in this model calculation. 

6 Summary 

In this paper the first measurements of the Balance Function in p+p, 
C+C and Si+Si interactions as well as centrality selected Pb+Pb collisions 
at ^snn = 17.2 GeV (the top SPS energy) are presented. 

The width of the BF decreases monotonically with increasing system 
size (from minimum bias p+p to central Pb+Pb collisions) by 24 ± 2% and 
with increasing centrality of Pb+Pb collisions (from peripheral to central 
collisions) by 17 ± 3%. A similar decrease, of the order of 14 + 2%, with 
centrality in Au+Au collisions was measured by STAR at ^Jsnn = 130 GeV. 
Thus the narrowing of the BF seems to be nearly energy independent from 
the top SPS to RHIC energies. 

Events from the string-hadronic HIJING model as well as shuffled events 
retaining only correlations from global charge conservation do not show any 
significant decrease of the BF width with increasing system size and cen- 
trality in nucleus-nucleus collisions. On the other hand, results from central 
Pb+Pb reactions at top SPS and Au+Au reactions at RHIC energies show a 
narrowing of the BF which suggests a delayed hadronization of the produced 
matter. For a more quantitative description of the data model calculations 
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have to include the effect of transverse flow of the matter at freeze-out. 

The energy dependence of the BF in the SPS range will be addressed in 
a future publication. 
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Interaction 


Number of events 


E Q range (GeV) 


(N w ) 


(b) [fm] 


P+P 


1M 




2 




C+C 


100K 




14 


1.9 


Si+Si 


100K 




37 


2.0 


Pb+Pb(6) 


300K 


29340 - 40000 


42 


11.5 


Pb+Pb(5) 


110K 


26080 - 29340 


88 


9.6 


Pb+Pb(4) 


88K 


21190 - 26080 


134 


8.3 


Pb+Pb(3) 


75K 


14670 - 21190 


204 


6.5 


Pb+Pb(2) 


100K 


9250 - 14670 


281 


4.6 


Pb+Pb(l) 


100K 


- 9250 


352 


2.4 



Table 1: Systems and centrality classes used in this analysis. Listed for p+p, 
C+C, Si+Si and six centralities of Pb+Pb collisions at 158^4 GeV are the 
range of the VCAL energy E , the mean number (Nw) of wounded nucleons 
and the mean value of the impact parameter. 



Interaction 


Year 


Pol. 


V x0 [cm] 


Ax [cm] 


Vyo [cm] 


Ay [cm] 


Vzofcm] 


Az[cm] 


p+p 


2000 


+ 


0.0 


1.0 


0.0 


1.0 


-580.0 


5.0 


C+C 


1998 


+ 


0.0 


1.0 


0.0 


1.0 


-579.1 


2.0 


Si+Si 


1998 


+ 


0.0 


0.3 


0.0 


0.5 


-579.5 


1.0 


Pb+Pb (m.b.) 


1996 




0.0 


0.1 


0.0 


0.1 


-578.9 


0.4 


Pb+Pb (m.b.) 


1996 


+ 


-0.05 


0.1 


0.05 


0.1 


-578.9 


0.4 


Pb+Pb (m.b.) 


2000 


+ 


0.0 


0.1 


0.0 


0.1 


-581.2 


0.4 


Pb+Pb (cen.) 


1996 


+ 


0.0 


0.1 


0.0 


0.1 


-578.9 


0.4 



Table 2: The different data sets used in the analysis. Listed for p+p, C+C, 
Si+Si and different sets of Pb+Pb collisions at y / J/v r /v = 17.2 GeV are the 
data taking period (Year), the field polarity (Pol.) as well as event selection 
cuts (see text for details). 
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y 


A [c/GeV] 


B [GeV/c] 


C [deg ■ (GeV/c) 1 / 6 } 


D [deg] 


-0.6 














-0.4 





-1 


63 


-8 


-0.2 








57 


-10 


0.0 





0.09 


63 


-13 


0.2 





0.08 


67 


-4 


0.4 


-7 


0.08 


65 


-3 


0.6 





0.05 


27 





0.8 








35 





1.0 





0.1 


41 





1.2 


0.34 


0.43 


109 





1.4 


0.36 


0.43 


100 





1.6 


0.55 


0.4 


100 





1.8 


0.6 


0.4 


88 





2.0 


0.61 


0.35 


73 





2.2 


0.73 


0.34 


55 





2.4 


1.7 


0.28 


60 





2.6 


2.8 


0.25 


60 





2.8 


5 


0.2 


57 





3.0 


7 


0.15 


60 





3.2 


7 


0.1 


70 






Table 3: Values of the parameters A, B, C and D of the acceptance curves 
(Eq. EJ). In the first column, the lower limit of the rapidity interval, y is 
given, y is calculated in the center of mass system assuming pion mass for 
all particles. 
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Interaction 


(Arf) (Data set 1) 


(A??) (Data set 2) 


(A77) (Data set 3) 








767 ± 007 


c+c 






721 ± 01 Fi 


Si+Si 






698 ± 01 1 


Pb+PbfVeto 6) 


698 ± 022 


695 ± 019 


704 ± 016 


Pb+Pb(Veto 5) 


0.695 ±0.022 


0.700 ±0.021 


0.689 ±0.021 


Pb+Pb(Veto 4) 


0.653 ±0.021 


0.672 ±0.019 


0.663 ±0.019 


Pb+Pb(Veto 3) 


0.642 ±0.021 


0.661 ±0.018 


0.645 ±0.019 


Pb+Pb(Veto 2) 


0.594 ±0.012 






Pb+Pb(Veto 1) 


0.582 ±0.011 







Table 4: The width of the BF for the three different data sets described in 
the text. 



Interaction 


(A77) (SHUFFLING) 


(A77) (HIJING) 


P±P 


0.784 ±0.007 


0.764 ±0.005 


C±C 


0.815 ±0.014 


0.746 ±0.010 


Si+Si 


0.833 ±0.011 


0.732 ±0.012 


Pb+Pb(Veto 6) 


0.823 ±0.020 


0.726 ±0.022 


Pb±Pb(Veto 5) 


0.823 ±0.021 


0.732 ±0.014 


Pb±Pb(Veto 4) 


0.806 ±0.021 


0.744 ±0.016 


Pb+Pb(Veto 3) 


0.804 ±0.022 


0.729 ±0.016 


Pb+Pb(Veto 2) 


0.807 ±0.015 


0.747 ±0.015 


Pb+Pb(Veto 1) 


0.818 ±0.018 


0.746 ±0.014 



Table 5: The width of the BF for the shuffled and HIJING data sets. 
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Figure 2: The acceptance curves in the pr — 4> plane for 2.5 < y < 2.7 (left 
plot) and 4.1 < y < 4.3 (right plot) at ^/snn = 17.2 GeV (color online). 




Figure 3: The pseudo-rapidity distribution of the accepted charged particles 
in central Pb+Pb collisions at y's/viv =17.2 GeV. 
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Figure 4: The BF versus A?7 for different centrality classes of Pb+Pb col- 
lisions for real data as well as for shuffled, mixed and HIJING events. The 
curves show Gaussian fits (color online). 
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Figure 5: The BF versus At] for real, shuffled and HIJING events (together 

with the Gaussian fits) for Si+Si (upper panel), C+C (middle panel) and 

p+p (lower panel) collisions (color online). 
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Figure 6: The width of the BF as a function of the event selections cut Az 
(upper panel) and track selection cuts | c?^ | marc (middle panel) and Iceman 
(lower panel) for p+p, Pb+Pb peripheral and Pb+Pb central data (color 
online) . 
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Figure 7: The BF versus At] for three centrality classes of Pb+Pb collisions 
together with the Gaussian fits. The distributions were normalized to the 
same integral for this comparison (color online). 
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Figure 8: The dependence of the BF 's width on the number of wounded 
nucleons for p+p, C+C, Si+Si and Pb+Pb collisions (color online). 
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Figure 9: The dependence of the BF 's width on the normalized impact 
parameter (b)/b max for Pb+Pb collisions (color online). 
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Figure 10: The dependence of the BF 's width on the normalized impact 
parameter (b)/b max , as measured by NA49 for Pb+Pb collisions at ^/snn = 
17.2 GeV and by STAR for Au+Au collisions at y/s^N = 130 GeV (color 
online) . 
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